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This article investigates the interaction of a spherically symmetric massless scalar field with a 
strong gravitational field. It focuses on the propagation of waves in regions outside any horizons. The 
two factors acting on the waves can be identified as a redshift and a backscattering. The influence of 
. . . backscattering on the intensity of the outgoing radiation is studied and rigorous quantitative upper 

bounds obtained. These show that the total flux may be decreased if the sources are placed in a 
region adjoining an apparent horizon. Backscattering can be neglected in the case 2mo/R « 1, 
that is when the emitter is located at a distance from a black hole much larger than the Schwarzschild 
radius. This backscattering may have noticeable astrophysical consequences. 
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I. INTRODUCTION. 



> , 

\ A black hole can be investigated by observing the gravitational radiation or the electromagnetic radiation emitted 
04 . by a medium surrounding it. The two major effects that act on the radiation are redshift and backscatter. While the 
redshift is well understood, little attention has been focused to date on backscattering, a phenomenon that prevents 
radiation from being transmitted exclusively along null cones. Backscattering can influence observations, since sharp 
and strong initial impulses can be weakened and dispersed. Electromagnetic or gravitational radiation, for instance, 
, reaching an asymptotic observer, can be weaker than expected. 
' In this paper we assess analytically this effect by considering as a toy model a spherically symmetric massless scalar 
field. This simplifying assumption should not restrict seriously the validity of our conclusions. We study the evolution 
I of pulses of radiation which are initially purely outgoing and emitted from a region close to a spherical black hole, 
^jpj Rigorous upper bounds on the magnitude of the backscattering effect are obtained; these estimates, while they break 
• • down close to the horizon, allow one to recognise the regions in which the backscattering may be neglected. We 
estimate how significant this effect would be in the case of neutron stars or black holes. The analytic estimates that 
are obtained here can be used in numerical gravity in order to assess numerical codes. They can be used to test both 
^ the accuracy and the long-time stability of the numerical models. The final section outlines possible astrophysical 
consequences of backscattering. 

Our results give the first, we believe, quantitative measure of this strong field effect. It is worthwhile pointing out 
that both our approach and the physical problem dealt with are different from that discussed in the existing literature 
on backscattering, which concentrates on the analysis of stationary radiation on a fixed geometric background 0| . 

II. SPHERICALLY SYMMETRIC MASSLESS SCALAR FIELDS IN MINKOWSKI SPACE 

Let us consider Minkowski space with a spherically symmetric coordinate system and a metric Qf^i, = 
(—1, +1, i?^, aiv? 9). Let t) be a spherical field which satisfies the simplest massless wave equation 

V''V^</) = 0. (2.1) 



The general solution (except at the origin) to Eq.(2.1) is 



f{R + t)+g{R-t) 
R 



(2.2) 



where / and g are essentially arbitrary functions of one variable. Therefore the general solution can be regarded as a 
superposition of an ingoing null field / and an outgoing null field g. Associated with the field (/) will be a stress-energy 
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tensor T^,y, an object of some importance, as we want to couple the scalar field to gravity.' A priori, it is required 
only that T^^ be a symmetric tensor that satisfies the conservation equation T^^ — 0. The standard choice is 



(2.3) 



This tensor will be subsequently used, but the reader should be aware that it is not necessarily unique. 

Given the way that the field runs along the ingoing and outgoing null directions, it is natural to consider the two 
ingoing and outgoing null derivatives of (f> 



When written in terms of U and V , the stress-energy tensor take a form 





^1{U^ + V^) =p 


Tqr 




Trr 




Tee 




T^<t> = 





J 



(2.4) 
(2.5) 



(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



From Eqs.(2.6) and (2.7) one can be tempted to interpret and as, respectively, outgoing and ingoing components 
of a radial flux. There is, however, a discrepancy with what can be' i nferre d fro m (. 2.2 ) ; from that follows, for instance, 
that field is purely outgoing if / = 0. When returning to Eqs.(2.4) and (2.£), it is clear that then U does not vanish 
identically. As a consequence, if the field is purely outgoing, the energy-momentum has an inward part which only 
vanishes asymptotically so that the energy-momentum is timelike even though the field is null. 

There do exist objects which more naturally reflect the splitting given by Eq.(2^) into ingoing and outgoing modes. 
Let us define 



RV 



h_{R,t)^^{dn + dt){Rcl,)=f' = ^ + ^. 
Extending R to cover the whole real line and defining h{R, t) via 



(2.11) 
(2.12) 



h{R,t) = h+{R,t),R> 
h{R,t) = h_{R,t),R< 



(2.13) 
(2.14) 



allows one to write the field equation (2.1) compactly as 

{dn + dt)h{R,t) = . 



(2.15) 



By comparison, when writing earlier the field equation in terms of U and V we get a pair of coupled equations. There 
exists another stress-energy tensor, more naturally related to h±, of the form 



T' 



T' — 



hi 



i?2 



i?2 

I 

Taa — T±± = 



Trr. — Tqq 



(2.16) 

(2.17) 

(2.18) 
(2.19) 



This also satisfies T^J^^ = in flat space. Unfortunately, it is not clear how to generalize that object to curved space 



so we will restrict our attention to the T^'^ as defined by Eqs.(|]^ - 2.10| ). While the field equation is best expressed 



in terms of h{R,t), as will be shown later, the conserved quantities are natural in terms of U and V. 
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The total energy inside a sphere of radius i? at a fixed time t reads 



r^{U^ + V^)dr 



It is now easy to compute various transport equations for the energy M 

dRM{R, t) = AttR^p = 7:R'^{U^ + V^) 
dtM{R,t) =-47ri?2j =ttR'^{U^ -V^) 
{dR + dt)M{R,t) =27ri?2[/2 
idt-dR)M{R,t) = -2T:R^V^ . 



(2.20) 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



It is clear from Eqs.(2.22) and (2.24) that and can be interpreted as the energy fluxes in the outgoing future 
and ingoing future null directions respectively. 



III. SPHERICALLY SYMMETRIC MASSLESS SCALAR FIELDS IN SCHWARZSCHILD SPACETIME 



Let us now consider a spherically symmetric massless scalar field minimally coupled to the Schwarzschild spacetime 
with the standard external metric 



ds^ = - 1 - 



2toc 
R 



(3.1) 



where mo is a positive constant, the Schwarzschild mass. While there is no explicit solution akin to Eq.(2.2), much 
of the analysis of Section 2 carries over. Define a 'radiation amplitude' h{s,t), se(— oo,oo), as follows 



h+iR,t) = h{R,t) =. h-ido + dR)iRcj)) 
2 7 

h^{R,t) = h{-R,t)= h-do + dR){Rcj)) . 
2 7 

where 7 = 1 — 2mo/i?. Given h{R) the field (j) can be recovered by a simple integration 



(3.2) 



1 



1 



^ = -^1 dr[h+{r)+h^ir)]^-^. 



(3.3) 



The scalar field equation, Vp9''0 = 0, Eq.(2.1), can be written as a single first order equation on a 'symmetrized' 
domain — oo < R < oo, 



It is useful to change coordinates and introduce the Regge- Wheeler coordinate [|j as a new radial variable 

r* =R + 2mo\n{- 1). 

2mo 



Using this, equation reads 



(3.4) 



(3.5) 



(3.6) 



From Eq.(|3.5D follows that dR/dr* — 1 ~ 2mQ/R and this allows us to write Eq.(3.6) as 



2moR 2mo - 

13-[(1 —'f^l 



R[ 



R 



(3.7) 



This can be solved by 
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i? 



(3.8) 



where the integration in Eq.(3.8) extends along the 'outgoing' null ray from a point on the initial slice at (r* — t,0) 
to (r*,t). (Of course in the negative R region, where we are dealing with the integration contour follows along 
the physical ingoing null ray.) It is clear that h m ixes the outgoing with the ingoing modes and thus is responsible for 
the backscattering. The leading term in Eq.( |3.8[ ) is just the standard gravitational redshift. One gets from Eg. (^?^) 
at t = 



ho{x) 



and ignoring the second term in Eq.(3.5) we arrive at 



h{r*iO) + T,T) = 1- 



R{x) 



1 - 



2mo 



i?(r*(0) +r) 



%*(0),0) , 



(3.9) 



(3.10) 



when r becomes large. 

As in the Minkowski space-time case the stress-energy Eq.(E^) naturally expressess via U and V defined below 



V = (7^5^ - 7-^^4)0 . 

The relevant components of the stress-energy tensor read 

The 'mass function' - the energy associated with a ball of radius R - is given by 



m{R,t) = An f r'^pd: 
Jo 



This is the analogue of the Misner-Sharp-Hawking-et al mass. Differentiation of Eq.(3.16) yields 



dnmiR, t) = AnR^p = TrR^{U^ + V^) 
dtm{R,t) = -47ri?2j= TT-fR^{U^ - V^) 



(3.11) 
(3.12) 



(3.13) 
(3.14) 
(3.15) 

(3.16) 



(3.17) 
(3.18) 



The derivatives of m along the outgoing and ingoing null directions can be written, employing Eqs.(3.17) and (3.18) 

as 



{-/dR + dt)m{R,t) -27r7i?2[/2 
{dt - "fdR)m{R,t) = -27r"fR^V^. 



(3.19) 
(3.20) 



Let us comment on the derivation of Eq.(3.1q), since it is not obvious and since that would help to understand why 



m{R,t) is chosen as above, in a non-canonical way. There is a timelike Killing vector = (1,0,0,0) which satisfies 
i(fi;v)- This implies, together with the stress-energy conservation law Tjf.^ = 0, (T'^C'^);^ — 0. The vectorial divergence 
can be written as {\/^-gTj^£^'^)^^ = 0; this leads to 



(3.21) 



The integral of the left hand side of Eq.(3.21) over a sphere is clearly dm/dt from Eq.( ^.16 ). The i ntegr al of the right 
hand side is just a flat space total divergence which becomes a surface integral. From this Eq.( 3.18 ) immediately 
follows. 
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IV. THE SPHERICALLY SYMMETRIC EINSTEIN - MASSLESS SCALAR FIELD SYSTEM 



In this section we will study the full Einstein-scalar field system, including the backreaction that the matter exerts 
on the geometry. The polar gauge condition trK = is our slicing condition. This guarantees that the normal to 
the slices is along the lines of constant R. Its major attraction is that to a large degree it eliminates the extrinsic 
curvature from the field equations. This allows one to express the metric explicitly in terms of the matter-field. 
The great disadvantage is that one cannot cope with horizons in the polar gauge, thus the latter is not suitable for 
the descriptions of late stages of the collapse but it is useful in the exterior zone that we wish to investigate. The 
aforementioned choice of gauge is equivalent to choosing a diagonal line element 

= -P{R, t)j{R, t)df + ^^^^dR^ + R^{r, t)dn^ , (4.1) 
j(R,t) 

where t is a time coordinate, i? is a radial coordinate that coincides with the areal radius and dQ^ = d9^ + sir? OdcjP' 
is the line element on the unit sphere, < < 27r and < < vr. At spatial infinity (of an asymptotically flat 
spacetime) the metric coefficient gim tends to -1-1 and gu tends to —1. Thus (3 and 7 go to -1-1 at infinity. We adopt 
the standard convention that Greek letters range from to 3 for spacetime objects while Latin indices range from 1 to 
3 for spatial objects. The Einstein equations G^^ = R^u — gfj_^R/2 = SttT^i^ can be written as a 1 -I- 3 system of initial 
constraints and evolution equations Q|. Let Et be a foliation of the spacetime by Cauchy hypersurfaces, with gij the 
intrinsic metric and Kij the extrinsic curvature. In the convention of Wald 2NKij — +dtgij and trK = g"^^ Kjj. 

As before, the scalar field is massless and spherically symmetric, with the stress-energy given by Eq.( ^.3| ). The 
matter energy density is p = — Tq and the matter current density is J = — Joij//3. As above, we define a 'radiation 
amplitude' h{s), se(— 00,00), as follows 

h+{R,t) ^h{R,t) ^]-{--do + dR){R^) (4.2) 
2 7 

h-{R,t)^h{-R,t)^ ].{-do + dR){R4>). (4.3) 



2^7 



One can show that 



where 



2i? 

can be expressed in the following useful form |^ 



/?(i?) = e-^^(/r+0^('-^)^ , (4.4) 



h^-l^ r dr[h+ (r) + (r)] ^h. (4.5) 



(3dr. (4.6) 



In spherically symmetric gravity one knows that the gravitational field is kinematical, any real dynamics is in the 
matter field. There fore t he m etric can be expressed in terms of the matter. In the polar gauge this can be done 



exphcitly (see Eqs. (}4.4| - [4.6[) ). 

The scalar field equation V^9^(/) = can be written as a single first order equation on a 'symmetrized' domain 
—00 < R < 00, 

(do + idR)h=ih-h)^^ . (4.7) 

This equation, together with definitions of h, h, (3, and 7, encodes the all information carried by the Einstein equations 
coupled to the scalar field. Notice that drh{r) = dr[h{r) + h{—r)] — drdr(j> — 0, where 4> = R4>j therefore 
initial data, ho, of compact support must satisfy the condition 



00 



drho = 0. (4.. 
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It has been shown shown, using relations between metric functions and their symmetry properties, that if (^^) holds 
true then drh{r, to) = in the existence interval of a solution. The local existence of a Cauchy solution for the 
above system and a global existence in an external region bounded from the interior by a null cone have also been 
shown ( i, i). 

An alternative representation of 7 that will be useful is 



7(i?) 



2too 2m,ext{R) 



\R\ 



R 



(4.9) 



where mg is the asymptotic mass and niext is a contribution to the asymptotic mass coming from the exterior of a 
sphere of a radius \R\, 



m,xt{R) - 4^ / -,{[Kr) - hf + [h(-r) - hf)dr 



Parallehng Eq.(|3l|) we can write 



(ao +75fl)me,t(i?) = -87r72(/i_ - hf 
(90 - 79fl)me.t(i?) = 87172 (/!+ - hf 



(4.10) 



(4.11) 
(4.12) 



Let us stress that many of the equatio ns wr itten down in this section take on a quite special form in the polar gauge. 
In particular, the object defined in Eq.(4.1C) is the Hawking mass (or mass function) which in a general gauge depends 
on the matter current as well as on the energy-density. In the polar gauge it simplifies, since the two null expansions 
are equal and the term which depends on the current drops out. Let us point out again that h{R), R > 0, represents 
outgoing radiation while h{—\R\) is an ingoing scalar wave. 

We can substitute Eq.(4.9) into Eq.(4.7) and rewrite the scalar field equation as 



{do + ldR)h ^{h-h) 



-2m{R)R(3 



(4.13) 



where m[R) = mo — mf,xt{R) is the Hawking mass at a ra dius R. This is similar to Eq.( p.7| ) but now backreaction is 
taken into account. Therefore our interpetation of Eq.(3.7) as giving a 'red-shift' due to the h term on the right-hand- 
side (determined by the local mass function, m(R), rather than the Schwarzschild mass, mo) and a 'backscattering' 
due to the h term sh ould remain valid in some appropriate limit. 

In going from Eq. (^^) to Eq.(3.7) we found what was essenti ally an integration function, (1 — 2mo/i?), which 
allowed us eli minate the term in h on the right-hand-side of Eq.(3_^). It turns out that this trick can be repeated 
with Eq.(|4.13|). 

Let us define 



In 



1 - 



2m(i?) 
R 



2m{r)dr 



r'^{l — 2m{r)/r) 



(4.14) 



where the integral is taken along an outgoing null ray. From Eq.(4.12) follows that the mass function monotonically 
increases along an outgoing null ray and also that m{R) — > msRo (the 'local Bondi mass' that can be assigned to 
a future null cone starting from Ro). Consider a null ray which starts at i?o; then m(i?o) < m(r) < msRo- This 
immediately gives 



2m{Ro) 



< 



2m(r) 



r2(l - 2m{Ro)/r) " r2(l - 2m{r)/r) 
Integration of this equation along the null cone yields 

2TO(i?o) 2m{R) 



< 



2mBRo 



r2(l - 2mBRo/r) 



1 



R 



> 1 - 



R 



> 1 - 



R 



(4.15) 



(4.16) 



Since m{Ro) — s- msRo as i?o 00 we conclude that rh{R) — s- niBRo- m(i?) is positive along an outgoing null cone 
hence it is clear that 1 — 2rh{R)/R, which is going to be our redshift factor, is monotonically increasing. Note that 
there is no guarantee that 1 — 2m{R)/R is monotonic. 



Using Eq( 4.14 ) and Eq.(19), it is easy to rewrite Eq.( 4.13 ) as 
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{do + -fdn){l 



2rh{R) 
R 



)h = h(l - 



or as a pair of equations on the half-line, i? > 0, 



R 
2m{R) 

R 



2m{R) 2m{R)R(i 



R 



2m(fl) N 2m(fl)/3 



(4.17) 

(4.18) 
(4.19) 



similar to equation (3.7). Reasoning as in Section III leads to a clean separation between redshift and backscattering. 



V. WEAK FIELD - REDSHIFT AND BACKSCATTER 

We can define the whole problem in the exterior domain, independent of whatever is happening in the interior. 
Already h, Eq.(|4.5|), is defined in a proper form. Given (4.8), h can be written as 



1 






2R 


J —OQ 





drh{r) 



(5.1) 



It is cle ar, f rom Eq.(4.4), that /3(i?) is only dependent on external quantities. The calculation of j{R), requires the 
use Eq.(4.6), together with the fact that, at infinity, 7 behaves like 1 — 2mo/R while /3 goes quickly to 1. This allows 
us to write 



7(i?) = 1 - 



1 



2r7io 



R R 



[1 - P{r)]dr ; 



thus there are two equivalent expressions for 'mg^,j.t{R), Eq.(4.10) and 

mext{R) = 



[1 - /3(r)]dr 



The form of Eq.( 4.10 ) inspires one to write Eq.( [4.4[ ) as 



I3{R) 



Eq.(4.9) gives the inequality 



7r r — 2m(r) ~ r — 2mo 



< 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



The factor 13 /^r can b e tak en out of the integral in Eq.(5^) and replaced with its value at r = i?. The remainder 
is then essentially Eq.( 4.ld| ). Select an e and an Ra such that simultaneously iti^Ra) ~ too, Tnext{RA) I mo < e and 
Ra> {1 + e)2mo. U R> Ra then 



1 > /3(^) > P{Ra) > e ~ 
In the same vein, using Eq.( |4.9|) we get 

7(i?) ^ 1 



"0 > e 



"0 



1 - Oimext/emo) 



2mo 

w 



+ O{mext/emo) 



(5.6) 



(5.7) 



Thus e and Ra should be chosen in such a way that mext/ewo <C 1. The last two expressions imply that the condition 
that backreaction is negligible is not that niext *C mo but rather mext ewp; in region s clo se to a horizon even if 
e ^ 1, a small cloud of matter can strongly influence the geometry. Estimates (|5.6|) and ( ^.7|) hold, for instance, for 
Ra > 2too(1 + \/mext/mo). In this case 



1 > KR) > mA) ^ 1 - 0{^m,xtlmo) , 



(5.8) 
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and 



7(i?) 



1 - 



2mn 



R 



0{y/ rriext/mo) ■ 



(5.9) 



In what follows our interest is focused on the region R > Stoqi since only in that region of spacetime we get sensible 
analytic estimates. This is equivalent to the c hoice e > 1/2. 

For R > 3mo the scalar wave equation Eq.(4.17) can be approximated as 



, o \fi In 2mo moR 

i9o + ldnKl~—)h^Hl-—)-— 



(5.10) 



from the preceding analysis it follows that the error terms are of order TOext/^o- Thus, in the limit of mext/'mo ^ 1 
we recover the Schwarzschild background scalar field equation, Eq.(3_^), with a solution Eq.(li) consisting both of 
the redshift and the backscattering terms. 

It turns out that mext can be used in another role. It is essentially the integral of the square of the first derivative 
of (f) and so one can use it to bound the pointwise value of = 2h. One can show, improving a coefficient in an 
inequality of ||], that 



1^1 < 



y/mext 



i?l/2 /87r(l - ^) 



(5.11) 



(An estimate similar to this, but written with the total Bondi mass mp in place of iTif.^t^ appears in [^.) Thus the 
solution of Eq.(5.10) can be estimated by solutions of following equations 



(5o + 9..)((] 



2mo 



)h 



2mo(l 



i?5/2j8^(l 



2mo 



(5.12) 



(ao + a..)((] 



2TOC 

R 



)h 



2mo(l- ^)V77w 
R^/^Js^il^^-^) 



(5.13) 



where the min us a nd plus signs gives a bound from below and above, respectively. The solution of Eq. ( ^.12 ) the 
same as of Eq.(3.8) 



_ 2mo \ 
R[r)) 



h{r*,t) = ho{r* -t) + 



Stt 



2mo 



and equation ( ^.13 ) is solved by 



Stt 



-2mo 



r,t) 
(r-t,Q) 



2mo 



2mo 



dv- 



1 



i?V-R- 2mo 



dv- 



1 



(5.14) 



(5.15) 



Notice that that the du in the integrals in Eqs.(5.14) and ( 5.15| ) is essentially dr* , the Regge- Wheeler coordinate 
introduced in Eq.(|3.6|) with dR/dr* = 1 — 2mQ/R. Therefore the integral can be written in closed form as 



rriext 



IT 



-mo 



dR 



R^y/{R-2mQ 



V_R-2mo I 1 

2mo-R 2mo\/2mo 



arctan 



R-2mo 



2mo 



167rrno 



2mo 
R 



1 - 



2mo ^ 



arctan 



_R-2mo 
2mo 



(5.16) 
(5.17) 



From Eq.( 5.17 ) not only can we see that the redshift term in the solution (Eq.( 5.14 ) or ( 5.15 )) is correct up to 
terms of order rriext /ino, we can also see that the total backscattering is bounded by a term of order niext/'m-o- 
Thus in the limit where mext/ino is small, the standard gravitational redshift represented by Eq.( 3.10 ) is recovered 
and this is valid along an outgoing null ray to an accuracy of Oy^uiext/'nT'O. Let us stress that this holds in a strong 
gravitational field. We need not assume R^ niQ, but only R > 3mo. 
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The estimate of the backscattering, Eq.(5.17), is amazingly regular. It equals when R = 2mo, remains always 
positive, monotonically increases and has a maximum value of i/Trmeait /64mo when R ^ oo. It rises quite rapidly 
to about 1/3 of its final value when R — 3mo and about 0.8 of its final value when R = 4mQ. This means that the 
backscattering is only important deep in the potential well and its contribution rapidly diminishes as one considers 
radiation which is starting off further and further out. The redshift factor, 7=1 — 2mo/R is much more severe and 
the energy flux (or the magnitude of h) is severely diminished if one starts close to the horizon. 

More detailed information about the nature and magnitude of the backscattering term will be given in the next 
section. 



VI. BACKSCATTERING 



Let us assume initial data which (at t = ) represent a pure outgoing wave, i. e., /i_ = or (equivalently) 
h{R,0) = for i? < 0. Let us take a point to the future {Ri,Ti), outside any horizon, Ri > 2mo; i?i > 3mo would 
even be better in order to guarantee that bacreaction is absent. Consider the ingoing future radial null ray which 
passes through (i?i,Ti). This will start from the initial hypersurface at some point (i?o,0) with _Ro > Ri- Along this 
null ray R monotonically decreases while mext monotonically increases. Eq.(5.17) implies 



Ri 



2mo / -I _ 2mo \ 
Ro I -Ro I 



ia M _ 2mo ) 
1 V J 



arctan 



arctan 



i?.o— 2mo 
2mo 



-Rj — 2mo 
2mo 



(6.1) 



where ruext is mext{Ri,Ti). This expression is valid to an accuracy of 0{'mext/n^o) and has an upper bound of 
a/ Trmei^t /64mo which is achieved when when Ti is large and i?i ~ 2too. It reduces rapidly as i?i increases, it is down 
to 1/5 of this value when Ri = 4mo as already discussed in Section 5. Far out along an outgoing null cone, where 
i?i ~ Ti ^ mo and where Rq w 2i?i, it goes to zero like R^ . In this regime the analysis of Sections 2 and 3 
would hold. Therefore one could expect to approximat e h by g{R — t) /2R. Further, approximating g{R — t) by a step 
function of finite width and substituting this into Eq.(3.8), one expects /i_ to fall off like R^^ . 



Take an outgoing null ray through the point (i?i,0) which forms the inner boundary of the outgoing wave. We 
estimate h- along this null ray, using Eq.(6.1), and h. The integration of Eq.(4.12) along this null ray yields an 



estimate of the total energy flux across this surface in the inward direction. This will be the total energy loss from 
the outgoing wave due to backscattering. 

To get explicit estimates, let us further approximate the integ ral in Eq.(^.l(S| ): R monotonically decreases along the 
ingoing lightray, so one can just replace the -\/l — 2mo/R by -y/l — 2mo/Ri. This yields 



mext 



-mo 



dR 



27r i?V(-R-2mo 



Thus we arrive at 



\lh-iR)\ < 



< 



< 



9 fdRR-^/^ 

27r(/ti— 2mo) ^ •f 



27r(flT-2mo)'^0 



3R, 



3R, 



3j2 



3/2 



27rmo 



3/2 



Ry/R - 2mo 



(6.2) 

(6.3) 
(6.4) 

(6.5) 



Knowing one can use its definition, Eq.(4.3), to compute Rtp and hence (/) and h. Eq.(4.3) yields the inequality 



(do + dr'){Rh) ^ -fh^ < 



27rmo 



TO, 



3/2 



i?(i?-2TOo)l/2 



(6.6) 



a similar inequality with minus sign gives a lower bound. 

These can be solved in the standard way to give (remembering to change from dr* to dR) 
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\Rh{R,t)\ < |i?/i(i?2,0)| + - 



ITlext 3/2 



dR 



27rmo "^ J (i?-2mo)3/2 
where we assume that the outgoing nuU cone starts at (i?2 , 0) . The integration gives 



\Rh{R,t)\ < \Rh{R2,0)\ 



3/2 



27rmo (i?2 - 2mo)i/2 



rUext 



3/2 



27rmo (i?- 2too)^/^ 



(6.7) 



(6.8) 



Now let us choose this outgoing null ray inside the outgoing burst of radiation; there is neither an ingoing nor an 
outgoing field, i.e., h+ — = 0. This does not imply (/)/2 — h = 0. Rather it only implies R(j) = C, with C a 
constant. This looks like the static field of a point charge. If there is a regular center, this would be excluded, but in 
the case of a black hole, or if there is some sort of complicated interior this cannot a priori be ruled out. On the other 
hand such a 'static' field has energy density and this will cause the mass function to vary along the lightcones. We 
are not interested in such variation focusing only in evaluating the change in the mass function due to backscattering. 
Therefore we make the extra assumpt ion that 0(i?2, 0) = 0. This means that the first ter m in Eq. (|6.8| ) vanishes. 
It is clear that the last term in Eq.(6.8) is strictly larger than \h- \ as given by Eq.(6.5) if i? > AniQ (i. e., when all 



radiation is placed outside 4too). Thus 

\h-h-\ < \h\ 



4 

< - 
- 3 



rriext 



3/2 n 

mJ R2 



27rTOo R{R2 - 2mo)i/2 ' 
The integration of Eq.(4.11) from Ra to 00 gives a bound for the total change in rriext-, 

16 f 2mo \^ f 1 — mo / R2 



Artip 



9 V R2 



1 - 2mo/i?2 



(6.9) 



(6.10) 



This expression demonstrates how sensitive the amount of backscattering is to the location of the innermost null 
cone. This estimate becomes meaningless if R2 ~ 3.5mo because Arriext ^ "mext- In the case of a neutron star, where 
2mo/i? < 0.1 (on the surface of the star) we have an upper bound for the backscattered energy of 2% of rriext- 

The above estimation can be significantly refined when initially outgoing pulses are far enough from the apparent 
horizon. The basic point uses the following argument. Take an outgoing null ray which starts at a point (-R2,0) and 
which goes through (i?i,Ti). Consider also the incoming null ray which starts at (Rq^O). These have corresponding 
tortoise coordinates r* (2), r* (1), r*(0) as defined by Eq.(p^). They satisfy 

r*(0) - t(1) = r*(l) = r*(2) + t(1) . (6.11) 

While it is NOT true that Rp < 2 i?i in general, it is almost true and it is certainly true, for instance when R2 > 12mQ. 
In going from Eq.(6^) to Eq.(6^) one can legitimately replace Rq by 2i?i. That allows one to get in (|6.5| ) the additional 
V2-1/2 



factor a 



V2 



(6.5) reads 



\Rh{R,t)\ < |i?/i(i?2,0)| + 2 



2^/2- 1 
3%/2 



ITlext 3/2 

mo' 

Trmo 



dR 



(i?-2mo)3/2 



(6.12) 



assuming that the outgoing null cone starts at (i?2, 0). While it is possible to solve this integral explicitly, an interesting 
approximation can be found by extracting the 1 — 2mo/i? and replacing it with 1 — 2mo/R2 to get 



\Rh{R,t)\ < \Rh{R2,0)\ 



8\/2-4 



rUext 



3/2 



3^2 V 2^mo (1 - 2mo/i?2)3/2 



1 



1 



Now the second term is more than twice as large as the estimate for \h^\. Thus we may conclude that 



\h-h\ < 



rUext 



3/2 



8V2-4 

3V2 V 27rmo (1 - 2mo/i?2)3/2 



1 



1 



RVIh 2R-2 



(6.13) 



(6.14) 



Therefore, the total change in rUext can be obtained by inserting the above into Eq.(4.11) and integration from R2 to 
00, which yields 



Q;^16/2mo 
Ame.,<^(— 



2(l^moAR2), 
1 - 2mo/i?2 



■rriext 



(6.15) 
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That estimate is better than the former one by the factor of o? ~ 0.4. Applying that to case of neutron stars, when 
the conditions assumed above hold true, we can get that the maximal amount of (possibly) backscattered radiation 
does not exceed 1 percent. 

The initially outgoing field, , generates a weaker ingoing field, which enters the 'no-radiation' zone behind 
the wavefront. This, in turn, scatters again off the gravitational field to generate a new outgoing field, which turns 
up at null infinity at a later time, after the first burst of outgoing radiation has gone. This is called the 'tail' term. 
Eg. (6. 9) give s an estimate for h in the 'no-radiation' zone. Let us substitute this back into the scalar wave equation, 
say Eg. ( 5. 10 ), and estimate the second order along the outgoing null ray. We get 



8 pm^ 



5/2 



3 V 27rmo \/R2 — 2too 



R r 



1 

i?3 



Vi?2 - 2too 



dR . 



Throwing away the second term in Eg. (3. If) and integrating Eg. (6. 16) gives 



\lh+{R,t)\ < - 



8 / mext 



5/2 



In the limit as R 



27rmo ^1 - 2mo/i?2 
oo along the null cone we get 



1 



1 



3i?: 



5/2 m^vm 



dR . 



|/i+(oo, oo)| < 



8 {mo/R2f/^ 



9 V 27rmo ^\ - 2mo/i?2 



(6.16) 



(6.17) 



(6.18) 



That should be compared with the leading term in ( |5.14 ), /iq; resorting to the definition of the external mass mext, 
we arrive to the conclusion that the tail term is smaller than the leading term by a factor (mo/i?o)^- 



VII. NUMERICAL RESULTS. 



In order to estimate guantitatively the effect of backscattering we solve numerically Egs.( 4.18 , 1.1£ and L5) ne- 
glecting backreaction (i.e. dropping t erms of order mext/'niQ) and calculate mext{R) using Eg.(4.1C) 



The basic evolution Egs.(4.18, 4.1£) are solved with a help of modified MacCormack predictor-corrector scheme 
|l6| . An integral defining h in Eg.( |4.5| ) is calculated by means of the extended trapezoidal rule. The whole procedure 
is second-order accurate in both space and time. 
We use the following form of initial data: 



h+{R)^Ao{R- 
h^{R) 



Ro)exp{-2Q{R- Rof/s^) 
= . 



(7.1) 
(7.2) 



whe re A q - amplitude, i?o ~ position of peak center, s - the measure of its width. This data fulfills the constraint of 
Eg. (HI). 

A sample of results is summarized in Figures 1 and 2 where the evolution of and h functions is shown. Notice 
the effects of the redshift and widening of the h+ profile caused by the light-cone structure in the vicinity of the 
horizon. 

As a guantitative measure of backscattering we take the mass flux ratio: 

(m 

ext 

{R, T) — m 

ext 



mextiRf.O) 



(7.3) 



where (i?j,0) and {R,T) are connected by an outgoing future radial null ray. 

Backscattering is noticeable close to the horizon in the regions where the h function is significant. We do not 
observe a significant effect outside an impulse - the observed flux through the inner boundary of a signal defined as 
a radius which contains initially 90 % of niext is less than 1 % (Figure 3, long dashed line). We do observe, however, 
a significant effect due to backscattering deep inside a region filled with radiation. The total leakage of mass is of 
the order of 10 % for Ri — Rq, where Rq is the position of the peak center in the initial configuration ansatz for a 
configuration starting from Rq = l.bRhor = 3mo (see Figure 4, solid line). This effect decreases with the distance 
from the horizon, for a similar configuration starting from Rq = SRhor = 6mo the effect is less then 2% (see Figure 3, 
short dashed line). 

From this we can conclude that backscattering deforms signals. Thus, assuming that a source mechanism for the 
radiation is known, we should in principle be able to detect backscattering in pulses or bursts of radiation coming to 
us from compact objects like massive neutron stars or black holes. 
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VIII. BACKSCATTERING: ASTROPHYSICAL IMPLICATIONS? 



We expect that massless spin-1 field (the photons) wiU behave in a very similar fashion to the massless spin zero field 
we have been discussing here. It is clear that two mechanisms effect the outgoing radiation. One, the redshift effect 
diminishes the intensity and frequency of the outgoing radiation, but the total energy in the radiation (as measured 
by the mass function) remains unchanged. This effect becomes larger and larger as the source approaches the horizon; 
the intensity is just multiplied essentially by 7 = 1 — 2mo/-R. 

The other effect is backscattering. The estimates that we derived in Section VI are not sharp but they are valid at 
least to within an order of magnitude, as the numerical results show. A quasi-stationary radiating plasma surrounding 
a black hole will produce radiation; some of the ingoing radiation will be reflected outwards (the black hole can act 
as a kind of a mirror) while some of the outgoing radiation will be reflected back. The reflection coefficients need 
not to be symmetric - presumably much more of the outgoing radiation will be reflected back - hence the radiation 
observed at infinity should be reduced below that expected when only the redshift factors are taken into account. The 
overall weakening of radiation in quasi-stationary systems may be of astrophysical significance in modelling activity 
of galactic nuclei powered by black holes, for instance a nuclei of the elliptical galaxy M87. That nuclei seems to be 
much less luminous than it should be, based on the standard accretion models ( ||^, ||^). Backreaction is certainly 
unable to to explain the full extent of the existing discrepancy, but it can contribute to the possible explanation if the 
bulk of the radiation is produced at distances of a few Schwarzschild radii, when backscattering can play an important 
role. 

Numerical results of Section VII demonstrate also another effect, namely up to 10% of the scalar field radiation 
emitted at i? = 3mo is lost from the main impulse - even if it reaches infinity, it does so with a significant delay. More 
precisely. Figure 3 shows, that for an extended pulse of radiation, its 'first' half is weakened, while the other gains 
in intensity. Thus the shape of the pulse is deformed. That opens another possibility for detecting backscattering, 
namely through the investigation of short-lived signals coming from an immediate vicinity of compact bodies. One 
class of phenomena, X-ray bursts, are already known to result from energetic processes on a surface of neutron stars 
PI , so they constitute a natural object of interest, albeit in that case the backscattering is rather weak, less than 
1 percent. The most interesting situation would be a merger of a neutron star and a black hole, since the effect 
would be then relatively strong. If gamma ray bursts |l^ result in the latter type of collision, then sufficiently precise 
observations should reveal imprints of bacscattering in the main burst as well in the afterglow. Their definitive absence 
would strongly favour alternative scenarios for the formation of gamma ray bursts. 
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FIG. 1. An example of h+ evolution. An initial profile defined by Eq.(7.f) with parameters Ao = 0.2, Ro = 3, s = 1, mo = 1 
(solid line) and profiles corresponding to i = 5 (long dashed), t = 10 (middle dashed) and t = 15 (short dashed) are plotted. 



13 



0.0000 



0.0005 



\ / 



/ 



-0.0010 

2^4 6 8 10 12 14 

FIG. 2. An example of h evolution. Parameters of initial confiuration and time sequence as in Figure 1. 
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FIG. 3. Examples of leakage through null cones - mass flux ratio rm{T) through different outgoing future radial null 
rays as a function of position along the ray R{T). The curves correspond to the following sets of parameters: solid line - 
Ri = Ro = 3,s = l,Ao = 0.2, mo = 1; short dashed line - Ri = Ro = 6, s = l,Ao — 0.2, mo — 1; long dashed line - 
Ri = 2.776, s = 1, Ao = 0.2, mo = 1. 
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